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ABSTRACT 

From the first publi s hed z > 3 survey of 21-cm absorption within the hosts of radio galaxies and 
quasars. iCurran et al.l (|2008bD found an apparent dearth of cool neutral gas at high redshift. From a 
detailed analysis of the photometry, each object is found to have a A = 1216 A continuum luminosity 
in excess of L1216 ~ 10 23 W Hz -1 , a critical value above which 21-cm has never been detected at any 
redshift. At these wavelengths, and below, hydrogen is excited above the ground state so that it cannot 
absorb in 21-cm. In order to apply the equation of photoionsation equilibrium, we demonstrate that 
this critical value also applies to the ionising (A < 912 A) radiation. We use this to show, for a variety 
of gas density distributions, that upon placing a quasar within a galaxy of gas there is always an 
ultra-violet luminosity above which all of the large-scale atomic gas is ionised. While in this state the 
hydrogen cannot be detected nor engage in star formation. Applying the mean ionising photon rate of 
all of the sources searched, we find, using canonical values for the gas density and recombination rate 
coefficient, that the observed critical luminosity gives a scale-length (3 kpc) similar that of the neutral 
hydrogen (Hi) in the Milky Way, a large spiral galaxy. Thus, this simple, yet physically motivated, 
model can explain the critical luminosity (L912 ~ £1216 ~ 10 23 W Hz -1 ), above which neutral gas is 
not detected. This indicates that the non-detection of 21-cm absorption is not due to the sensitivity 
limits of current radio telescopes, but rather that the lines-of-sight to the quasars, and probably the 
bulk of the host galaxies, are devoid of neutral gas. 

Subject headings: galaxies: active — galaxies: ISM — radio lines: galaxies — ultra violet: galaxies — 
galaxies: high redshift — cosmology: early universe 



1. INTRODUCTION 

Hydrogen gas accounts for 75% of all the baryonic 
matter in the Universe, of which the cool component, 
the reservoir for star formation, is traced by the radio- 
band 21-cm spin-flip transition. Due to the low prob- 
ability of the transition, compounded by the inverse 
squa re law, this is essent ially undetectable at z > 0.2 
(see iCatinella et al.l 120081 ). although in absorption the 
line strength is dependent only upon the column density 
of the absorbing gas and the radio flux of the background 
source. 

Hydrogen has been detected in the ultra-violet band 
Lyman-a transition, which traces all of the neutral gas, 
in 1500 high redshift galaxies intervening the sight- 
lines to more distarrt_2uasi^teUar_^ i iects (QSOs, sec 
ICurran et all 120021: iNoterdaeme et al.l l2009h . However, 
despite four decades of searches, knowledge of the cool 
component of this gas in the distant (z > 0.1) Universe 
remains very scarce, with only 42 cases reported in these 
absorbers, intervening radio-loud QSOs ( quasars j]], in 
addition to 35 associated with the quasar host galaxy 
itselffl 



In both cases, the majority of detections occur at red- 
shifts of z < 1 (look-back times < 7.7 Gyr)0 In the 
case of the intervening absorbing galaxies, the apparent 
lack of cold gas at high redshift may be accounted for 
by geometry effects: In an expanding Universe absorbers 
at redshifts of z It 1 are always disadvantaged, in com- 
parison to the low redshift (z < 1) absorbing galaxies, 
in how effectively the abs orber can cover the higher red- 
shift background source (|Curran fc Webbl[2006TlCurranl 
l20Tl . 

Since for the associated systems the absorbing gas is 
located within the quasar host galaxy, such geometry ef- 
fects cannot account for the fact that the 21-cm detection 
rate at z < 1 is double that at z > 1. Furthermore, only 
one as sociated 21-cm a bsorber has ever been found at 
z > 3 (jUson et al]|199li ). This runs contrary to the ex- 
pectation that at these redshifts (look-back times > 11.5 
Gyr), much of the gas has yet to be consumed by star 
formation, meaning that we would expect the abundance 
of hydrogen to be many times higher t han in the present 
day Universe (e.g. iPeroux et al.ll200H ). 

In addition to these covering factor effects, for a given 
column density, the optical depth of the 21-cm absorp- 
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1 Compiled in Curran (201 0|), w it h the addition of thos e recently 
reported by Srianand et al. (2010); Curran et al. (2011c). 

2 Compiled in Curran & Whiting (201(J), with the addition 
of three new associated absorbers, two reported in ICurran et al] 



|2011a |) and one in ICurran et al.l H2011dl) . See also lAllison et al.l 

3 We employ a standard A cosmology with Ho = 
71 km s" 1 Mpc" 1 , flatter = 0.27 and fi A = 0.73. 
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tion is dependent upon t he spin temperature of the gas 
(I Wolfe fc Burbided 119751) . Since only atoms populating 
the lower hyperfine level can absorb in 21-cm, the spin 
temperature may be elevated through: 

1. Excitation to the upper hyperfine level 
(jPurcell fc Fieldlfl956l) . 

2. Excitation above the ground state, particularly by 
Lyman-a (A = 1216 A) photons (|Fieldlll959f ). 

3. Ionisation. 

Although, with the data available, excit ation to the up- 
per h yperfine le vel through c ollisions ([Purcell fc Field! 
19561 see al s o iDraind I2011D cannot be ruled out, 



Curran et al.l (|2008bf) find no dependence of the 21-cm 



detection rate on the rest-frame 1420 MHz continuum 
luminosity of the active galactic nucleus (AGN), thus at 
le ast ruling out e xcitati on through this radiative process. 

iCurran et al.l (|2008b[) do, however, find a strong de- 
pendence on the rest-frame A = 1216 A (y = 2.47 x 10 15 
Hz) continuum luminosity. Specifically, that 21-cm ab- 
sorption has never been detected above an apparent crit- 
ical luminosity of ii 2 ie ~ 10 23 W Hz -1 . For a » 50% 
detect ion rate at L1216 

< 10 23 W Hz" 1 jC urran et al.l 
l2008bf ) . the probability of detections out of 19 searches 
occuring by chance is 1.9 x 10~ 6 (significant at 4.76g 
assuming Gaussian statistics, ICurran et al.l l20lTal) . So 
although the gas may be excited through other processes 
(collisions and the Cosmic Microwave Background), this 
correlation strongly suggests that excitation above the 
ground state (and possible ionisation) by A < 1216 A 
photons is the dominant cause of the non-detections. 

Given that 17 of the 19 L 12W > 10 23 W Hz" 1 sources 
are type-1 AGN, it is also possible that the absorption, 
by cool gas in the circumnuclear obscuring torus invoked 
by unified schemes, simply does not occu r along our line- 
of-sight to the continuum source (e.g. [M organt i et al.l 
[200lt iPihlstrom et al.l [2001 IGupta et all I2006D . How^ 
ever, at £1216 < 10 23 W Hz -1 , bo th type-1 and type- 2 
AGN exhibit a 50% detection rate (jCurran et al.ll2008bl 
I2011al) . indicating the absorption must primarily arise 
in the main galactic disk, which is randomly oriented 
with respect to the torus. Therefore, the bias towards 
type-1 objects at Lme ~ 10 23 W Hz -1 is due to these 
tending to arise in the more luminous quasars, as op- 
posed to radio galaxies, which tend to be associated with 
type- 2 objects and is the refore not an orientation effect 
(ICurran fc Whitindl2010l) . 

Thus. lCurran et al.l (|2008bh interpreted their exclusive 



non-detections at z > 3 to the traditional optical selec- 
tion of targets, in conjunction with the high redshifts, 
introducing a bias towards the sources more luminous 
in rest-frame ultra-violetQ The same critical ultra-violet 
continuum lumi nosity is also evide nt in the lower red- 
shift surveys (see [Mi ison et aLll2012f) and attributing the 
lack of cold gas in the hosts of powerful AGN to the high 
ultra-violet luminosities exciting the gas beyond detec- 
tion (which we dub "the UV interpretation' '), can explain 

4 Despite shortlisting the faintest objects (with blue magnitudes 
of B > 19, see figure 5 of ICurran et~aTll2u09T) . 



why this effect is seen at all redshifts. The UV interpre- 
tation may also account for several other issues in extra- 
galactic radio astronomy, such as the elevated detection 
rate in compact objects and the prefe rence for 21-cm de- 
tectio n in radio galaxies over quasars ([Curran fc Whiting! 
[20101) . 

Given the low probability of zero detections occuring 
by chance above a given 1216 A continuum luminosity, 
in conjunction the fact that A < 1216 A photons excite 
(and possibly ionise) the gas so that it cannot absorb in 
21-cm, there is little doubt that the UV interpretation is 
the correct physical description. This has been confirmed 
by an independent survey for 21-cm in 143 radio sources 
at redshifts 0.02 < z < 3.8, where the lack of detections 
is cor related with t he UV luminosity ([Grasha fc Darling] 
120111) . as well as bv lPage et al.l ([2012D . who find a critical 
X-ray luminosity, above which sources are not detected 
in 250 (im continuum emission, a tracer of star formation. 

However, one question remains unanswered: Why is 
there a hard limit to the UV luminosity, above which 
the gas is excited beyond detection by the most sensi- 
tive radio telescopes, rather than a continuum where the 
detections gradually become fewer and fewer as the ultra- 
violet luminosity increases? We address this issue here. 

2. PHOTOIONSATION EQUILIBRIUM 

For a cloud of hydrogen containing an ionising source, 
the equilibrium between photoionsation and recombina- 
tio n of protons and e lectrons in a nebula can be written 
as (|Osterbrocklll989j ). 



-7— dv = 47T 

hu 



n p n e aA,B r 



2 dr, 



(1) 



where L v is the specific luminosity at frequency v and 
h is the Planck constant, giving the number of ionising 
photons per second. On the right hand side, r; on is the 
extent of the ionisation, n p and n G are the proton and 
electron densities, respectively, and ola,b the radiative 
recombination rate coefficient of hydrogen (see Sect. l2.2[ ). 

Since, after excitation to the upper hyperfine level, the 
next excitation is to n = 2 by Lyman-a photons, our 
proxy has been the A = 1216 A continuum luminos- 
ity. However, since excitation to the n = 2 level and 
ionisation of the hydrogen atom are so close in energy 
(both events being « 2 x 10 6 times as energetic as the 
spin-flip transition), this critical luminosity should also 
apply in the case of ionisation. In order to verify this, 
in Fig. [1] we show the A = 912 A continuum luminosity 
distribution. The luminosities h ave been derived from 
the photometries as described in ICurran et al.l (|2008b[) . 
but with the inclusion o f data from the Ga laxy Evolu- 
tion Explorer (GALEX. [Martin et al.ll200l . These, in 
conjunction with the BVRK magnitudes from the liter- 
ature, allow reliable power-law fits to the rest-frame UV 
data (corrected for Galactic extinction using the maps of 
ISchlegel et al1li~998l) over a range of redshifts, from which 
the A = 912 A cont inuum luminosities were derived (see 
ICurran et aT1l20T2T) . 

From Fig. [TJ we see that the same approximate crit- 
ical value applies in the case of ionising photons. That 
is, 21-cm absorption has never been detected above a 
luminosity close to L912 ~ 10 23 W Hz^ 1 . The largest 
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Fig. 1. — The A = 912 A continuum luminosity (where available, see Curran ct al. 2012T) versus redshift for the sources searched in 
21-cm absorption. The filled symbols/hatched histogram represent the 21-cm detections and the unfilled symbols/unfilled histogram the 
non-detections. The shapes represent the AGN classifications, with triangles representing typc-1 objects and squares type-2s (+ and x 
designate an undetermined AGN type for a detection and non-detection, respectively). 
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measured 912 A luminosity for which there is a detec- 
tion is L912 = 1.1 x 10 23 W Hz -1 , above which there 
are 20 non-detections. Of the sources for which we could 
reliably determine £912, there are 38 detections and 60 
non-detections (i.e. a 39% detection rate) below this lu- 
minosity. Applying this probability of p = 0.61 for a 
non-detection to the L912 > 1.1 x 10 23 W Hz -1 sources, 
there is a binomial probability of 5.09 x 10 -5 of the 20 
non-detections occuring by chance, a 4.05<r significance. 

Thus, although only excitation above the ground state 
is required to explain the dearth of 21-cm absorption in 
UV luminous sources, it is possible that ionisation of the 
gas is the primary cause of the non-detections. Given 
that the lifetime in the n = 2 state is only <~ 10 -8 sec, 
this is the more likely situation and so we are justified in 
applying Eq. [T]to this problem. Thus, in Sect. 12.11 we 
derive the value of the left hand side of the equation, for 
the sources searched in 21-cm, and in Sect. 12.21 we apply 
various recombination models to the right hand side of 
the equation. 

2.1. Photoionsation rates 

In addition to determining the UV fluxes from the 
GALEX photometries and the BVRK magnitudes, in 
order to investigate differences between the 21-cm de- 
tected and UV luminous non-detected samples, we ob- 
tained multi-wavelength data from all of the relevant 
photometries given by the NASA/IPAC Extragalactic 
Database (NED), aga in correcting for Gal actic extinc- 
tion using the maps of ISchlegel et al.l (|1998f ). As per the 
A = 912 A continuum luminosities, above, we used the 
rcdshifts of the targets to blue-shift the observed frequen- 
cies back to the source rest-frame values and converted 
the observed fluxes to luminosities. We then averaged all 
of the luminosities within a specified frequency range to 
obtain a composite SED for all of redshifted the sources 
searched in 21-cm absorption. 

Performing a Kolmogorov-Smirnov test between the 
binned luminosities of the various sub-samples (shown in 
each panel of Fig. [2]), we find no evidence that the 21-cm 
detected and Lyv ~ 10 23 W Hz -1 non-detected samples 
are drawn from different populations, with a probability 
of > 0.2 (for all bin widths) that they are drawn from 
the same population. However, between either the 21- 
cm detected/Luv ~ 10 23 W Hz -1 non-detected samples 
and the Luv ~ 10 23 W Hz -1 sample, it is seen that, 
for sufficiently high resolution bins, the probability can 
get as low as 1 x 10 -6 , although this is due to the ex- 
tra high frequency points in the Luv ~ 10 23 W Hz -1 
non-detected sample, with a probability of « 0.004 be- 
ing more likely. This still suggests, however, that the UV 
luminous sources are drawn from a different sample than 
those with L UV £ 10 23 W Hz" 1 . 

In order to obtain the rate of ionising photons, we are 
interested in frequencies above v = 3.29 x 10 15 Hz. How- 
ever, as seen from Fig. [2 there is a large gap in the SEDs 
between ~ 10 16 and ~ 10 17 Hz, the range of spaced-based 
ultra-violet observations between the optical and X-ray 
bands. Although the X-ray observations are also space- 
based, these generally have more sky coverage than the 
ultra-violet observations and are thus more likely to have 
observed one of our sources. So, in order to obtain an es- 
timate of j v , (L v /v) dv, where v mVL = 3.29 x 10 15 Hz, we 



smooth the SEDs (Fig. [3]) and interpolate a power-law 
fit between ~ 10 16 and ~ 10 20 Hz to obtain the mean 
dependence of L„ on v. The photon rate is given by 

/ — - dv, where log 10 L„ = odog 10 v+C => L u = 10 c u a 

Ju. 

for a power-law, where a is the spectral index and C the 
intercept. Solving this, 

10 c r°° 10 c -10 c 

/ u a ~ l dv=^- [1/*]°° = -^-v a where a < 0. 
h J v ah ion ah 

From the composite SEDs, for the non-UV luminous 
sample we find L„ w 10 37 3 1/~°- 95 W Hz -1 , giving 
5.5x 10 55 ionising photons sec -1 and for the UV luminous 
sample, L u w 10 346 v~ Qm W Hz" 1 , giving 2.9 x 10 57 ion- 
ising photons sec -1 . This is f=a 50 times the luminosity 
of the 21-cm detected sample, which is consistent with 
the factor of « 7 in the luminosity distances between 
the z It 3 sample and the cluster of 21-cm detections at 
z < 0.9 (Fig. 0J. 

2.1.1. The critical photoionsation rate 

Since we are interested in the ionising photon rate re- 
sulting from a critical luminosity of Luv ~ 10 23 W Hz -1 
(Fig. [T]), we use the highest 21-cm detected luminosity 
of L912 = 1.1 x 10 23 W Hz" 1 at 3.29 x 10 15 Hz, together 
with the above spectral index of a — —0.68, to obtain 
L v w 10 33 - 6 is"" 68 W Hz -1 , which gives 2.9 x 10 56 sec" 1 
for the critical ionising photon rate. Referring to the 
literature, f r om th e spectra of several hundred QSOs, 
iTelfer et al.1 (I2002T) find a mean optical-X-ray slope of 
a = —1.5. This is significantly steeper than the mean 
spectral index derived for our sample, which consists ex- 
clusively of powerful radio sources, although applying a 
critical luminosity of L912 = 1.1 X 10 23 W Hz -1 gives 
L„ w lO 46 - 3 ^- 1 - 5 WHz- 1 => 1.1 xlO 56 sec" 1 , which is in 
the ballpark of the value derived for our sample. This is a 
consequence of the steeper spectral index being compen- 
sated by a larger constant (intercept) and the fact that 
the lower frequency end of the UV SED [y ~ 3.3 x 10 15 
Hz) contains most of the energy. We are therefore con- 
fident in applying J u {L v /hv) dv = 2.9 x 10 56 photons 

sec -1 to the left hand side of Eq. [TFI 

Lastly, it is clear that, while the specific continuum lu- 
minosity (L1216 or L912) may provide an indicator of the 
amount ionising radiation from the AGN, the integrated 
ionising luminosity (i.e. the ionising photon rate) is the 
correct measure. By fitting poly nomials to the phot om- 
etry of the individual sources (ICurran et al.l 12011) . 111 
Fig. @] we re-plot Fig. [T] in terms of J {L v /hv) dv. For 
those for which accurate polynomial fits could be ob- 
tained, we see that above ~ 10 56 ionising photons sec -1 , 
21-cm searches have resulted in exclusive non-detections. 
The highest photon rate for a detection, which could be 
reliably determined, is 1.7 x 10 55 sec -1 , a value above 
which there are 29 non-detections. Below this rate there 

5 Other studies of the UV continuum slope in quasars and AGN 
indicate their spectral indices to be in the range a = —1.4 to —0.56, 
although with considerable scatter. These give cr itical photon rates 
of ~ 5 x 10 55 an d ~ 2 x 10 54 - 2 x 10 57 sec" 1 flScott et al.ll200l 
IShuII et aI.l[20Ta. respectively). 
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Fig. 2. — The composite SEDs for each of the various sub-samples for a bin width of ±0.1 (i.e. each bin is 10 2 



1.6 times the 



frequency of the previous). In the left panel the vertical dotted line shows the source- frame frequency of v = 3.29 X 10 15 Hz (A = 912 A), 
above which the photons ionise the atom. The right panel shows the probability, from a Kolmogorov-Smirnov test, that the two samples 
tested are from the same population for various bin widths for both the whole frequency range (filled circles) and A < 912 A (hollow stars). 



6 



Curran & Whiting 



are 38 detections and 51 non-detections, giving ap = 0.57 
probability of a non-detection. Using this proxy, the 
binomial probability of 29 out of 29 non-detections oc- 
curing by chance is just 8.32 x 10 -8 . This 5.36c result 
therefore strongly suggests that ionisation of the gas by 
A < 912 A photons from the AGN is responsible for 
the non-detection of 21-cm absorption in high redshift 
sources. 

2.2. Recombination models 

We now parametrise the right hand side of the pho- 
toionsation equilibrium expression (Eq. [1} . Since we are 
concerned with the ionisation of neutral gas and its sub- 



sequent recombination, 



Also, in optical 



band observations of an optically thick plasma, where 
direct capture onto the ground state is excluded, qb is 
used. However, since we are concerned with the ground 
state, oa is the relevant total rec ombination rate coeffi- 
cient ([Osterbrock k, Ferl and 200(|) . We choose this value 
at 2000 K, a A = 1.27 x 10" 12 cm 3 sec"^ the typical 
upper limit to the spin temperature found in intervening 
absorbers (when the Lyman-a line is also detected and an 
upper limit to the spin temperature can be determined, 
iCurranet ani20To( ). 

Naively assuming a constant particle density of n = 
10 cm -3 (typical of the cool neutral 21-cm absorbing 
interstellar medium) throughout the nebula, we find 
Hon = 3 and 13 kpc for the mean radii of the "Stromgren 
spheres" of the UV non-luminous and UV luminous sam- 
ples, respectively. Although the latter value is of the 
same order of magnitude as the extent of neutral gas in 
a large galaxy, this model represents a gradual increase 
in ionised radius with luminosity, with no critical value. 

2.2.1. Exponential gas density distribution 

Unlike the idealised ionised region around a star, we do 
not expect the gas density to remain constant on galactic 
scales. A more realistic model of the density of the cold 
neutral medium (CNM) within a galaxy is that of an ex- 
ponential de crease in the gas density with distance from 
the nucleus (jBegeman et alj|1991l: Kalberla et al.ll2007t ). 
Thus, for n = n e~ r / R , where n is the gas density at 
r = and R is a scale-length describing the rate of decay 
of this with radius, Eq. Q] becomes 



— - dv = An a a n 2 n 
hv 



- 2r / R r 2 dr 



7r a a tIq R 3 — R e 



-2n 



»/-B(2i 



+ 2 r ion R + R 2 



(2) 



Unlike the constant density distribution, this be- 
comes independent of r at sufficiently large radii, i.e. 
J (L v jhv) dv — > 7r aA iip R 3 . Conversely, for a given 
scale-length, R, there always exists a "ceiling luminos- 
ity" (the number of ionising photons x h) for which all 
of gas is ionised (Fig. [5]). 

Using the above values of a a = 1-27 x 10~ 12 cm 3 
sec" 1 and uq = 10 cm -3 , the critical ionising photon 
rate of 2.9 x 10 56 sec -1 gives a scale-length of R = 2.9 
kpc. We can compare this to the H I in the Milky Way, 

6 Compared with ere 
of oa or «b making little difference. 



where iKalberla fc Kerd (|2009f ) fit an exponential pro- 
file to the mid-plane volume density distribution to find 
R = 3.15 kpc and n Q = Q.9e R °/ R = 13.4 cur 3 . This is 
in close agreement with our values, demonstrating that 
the mean SED normalised by a A = 912 A continuum 
luminosity of 10 23 W Hz -1 is sufficient to ionise all of 
atomic gas in a large spiral galaxy, rendering it unde- 
tectable in 21-cm. 

It is therefore clear that an exponential decrease in 
gas density with distance from the nucleus can naturally 
yield a critical value in the UV luminosity which is close 
to that found obscrvationally. This does however, rely on 
a simple model of the CNM, within which various struc- 
tures and phases will be embedded, such as the warm 
neutral medium, as well as localised regions of ionised 
gas and dense molecular clouds. However, here we are 
modelling the large-scale CNM, for whic h an exponential 
density distribution is a re alistic model (jBegeman et al.l 
fl99lt IKalberla et al .1120071) . 

Although using the canonical values for a and no gives 
the correct scale-length for the observed photon rate, a 
further physical (sanity) check can be obtained by deriv- 
ing the total jgas mass from the gas density and volume 
via M gas = J Q pdV . In this case, where the particle den- 
sity of n protons cm -3 corresponds to p = 1.67x 10 -21 x?i 
kg m~ 3 , exponentially decaying with r across a disk of 
thickness, t (r) = r//pL, we have 



27rno / er r/R rtdr 
Jo 

2irn R 3 ~ 



/fl 



2-7T uq 
/fl 

,2 



e -r/R r 2 







1 R 2 R 



,(3) 



where the flare factor, /fl, describes the flaring of the 
H I gas scale- height with galactocentric ra dius. Applying 
the mean Milky Way value of /fl ~ 20 (jKalberla et al.1 
120071 ). a scale-length of R = 2.9 kpc gives a total gas mass 
of M gas = 7.5 x 10 9 M Q (Fig.EJ bottom panel). This is 
close to the mean value found from a low redshift survey 
of 21-cm emission fro m the 1000 Hi bright est galaxies 
in the southern sky (|Koribalski et al.l 120041) , giving us 
further confidence in the exponential decay model and 
choice of gas density. 

2.2.2. Alternative temperatures and disk profiles 

Although spin temperatures in intervening 21-cm ab- 
sorbers may be, on average, ^ 2000 K, without the to- 
tal neutral hydrogen column density from the Lyman-a 
transition, as is the case for the associated absorbers, 
an upper limit to the spin temperature cannot be com- 
puted. Traditionally in the discussion of ga s ionisation, 
temperatures of ~ 10 4 K are brandished (e.g. lOsterbrockl 
Il989t lHaiman fc Reesll200l . although these are in the 
case of Lyman-a emission, rather than the much less en- 
ergetic 21-cm absorption. Given that the gas is most 
likely ionised however, in Fig. |6] (top panel) we show the 
ionising luminosity versus the extent of the ionised gas 
for aA = 4.19 x 10~ 13 cm 3 sec -1 (i.e. at 10 4 K ). From 
the observed critical rate of 2.9 x 10 56 photons sec -1 , we 
see that the scale-length increases to R — 4.2 kpc, which 
is larger than that of the Milky Way. Thus, for this tem- 
perature the photon rate is more than that required to 
ionise all of the neutral gas, while demonstrating that 
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our main result is not overly sensitive to the choice of 
temperature Q 

Although having no effect on the extent of the ion- 
isation, in order to investigate the effect of a different 
disk profile on the mass, in Fig. [6] we show the dis- 
tribution of gas mass for a sphere, rather than a disk 
(i.e. dV = Aitr 2 dr). For R = 4.2 kpc, this gives 
a total gas mass of M gas = 4.6 x 10 11 M©, which is 
close to the total (dynamical) mass expected in a galaxy 
and is thus too high@ and, applying a temperature of 
2000 K (i.e. R = 2.9 kpc), lowers this only slightly 
to M gas = 1.5 x 10 11 M Q . This confirms that the disk 
model, which reproduces a gas mass close to the typically 
observed value, is the more physically accurate distribu- 
tion. Furthermore, although the most luminous qu asars 
may reside in elliptical galaxies (jTavlor et al.lll996() , this 
shape traces the stellar distribution and not necessar- 
ily that of the neutral gas. For instance, "superdisks" 
of gas and dust in the elliptical hosts of powerful ra- 
dio galaxies have been prop osed (jAthreva et al.l 119981 : 
iGopal-Krishna fc Wiital 12000ft . wi th diameters of > 75 
kpc ( Gopal-Krishna &: Wiitall2000ft . perhaps up to w 300 
kpc (|Curran et al.ll2011dl) . 

2.2.3. Alternative gas distributions 

For completeness, we investigate several alternative 
density distributions for the gas. These are typically 
profiles that arise from dynamical models, often applied 
to the dark matter content of a galaxy, although we 
are interested in their effects when applied to the dis- 
tribution of the CNM. The Jaffe profile (fjaifel [l98l 
models the distribution of light in a spherical galaxy as 
n = no (r s /r) 2 /An (1 + r/r s ) 2 , where r s is the radius 
which contains half the total emitted light. In this case 
Eq. [1] becomes 



— - dv = 47T Of A «n 

hv 



(r s /r) 4 



; tt a A nl r 3 



(4tt) 2 (1 
r s (3 r s 3 + 22 r 2 r + 30 r s 



dr 



r/r s ) 4 
2 + 12r 3 ) 



r (r s 



r )3 



-121n(r s +r) + 121n(r s ) 



■(4) 



However, due to the "cuspy" nature of the distribution, 
solving this over these limits yields infinities and between 
any reasonable limits yields unreasonably large numbers, 
even when the approximation n = hq (r s /r) 2 /4 ir (1 + 
r/r s ) 2 « n (r s /r) 2 / An (r /r s ) 2 = (n /4Tr)(r s /r) 4 is used. 

A similarly asymptotic density distribut ion is given by 
the N avarro-Frenk- White (NFW) profile dNavarro et al.l 
11996ft . which models the density as the distribution of 
dark matter in the halo, via n = hq (r c /r)/(l + r/r c ) 2 , 
where no and r c are the core density and radius of the 
halo, respectively. Here the right hand side of Eq. [T] 
becomes 



4-7T qa 



2 2 

n r 



47T qa n 



(r c /r) s 



(l + r/r c )' 



■ dr 



-7ra A n 2 Q rl 



1 



(5) 



Unlike the Jaffe profile, the photoionsation equilibrium 
equation can be solved, again giving a ceiling luminosity, 
albeit less pronounced than for the exponential gas dis- 
tribution (Fig. top panel). The L uv = 10 23 W Hz" 1 
threshold of 2.9 x 10 56 A < 912 A photons per second 
gives r c = 2.6 kpc, which is very close to the scale-length 
of the expon ential disk and typic a l of that found f or near- 
by galaxies (|de Blok et al.|[200l IQh et al.|[20TlT) . For a 
spherical mass distribution the NFW profile gives 



M gas = 4?r n 



(1 



47T no ?"c 



- 1 + ln(r c 



r/r c ) 2 
r) + ln(r< 



' dr 



(0) 



from which r c = 2.6 kpc gives a total gas mass of 
Mgas = 7.1 x 10 11 M Q , which, not surprisingly given 
the distribution used, is close to the expected value for 
the dynamical mass. 

The asymptotic density distribution of the NFW pro- 
file can be avoided by employing the halo density dis - 
tribution of an isothermal sphere (jBegeman et al.l[l991[ ) , 
i.e. n = ?io/(l + r / r c) 2 , Fig. [8] (middle panel). For this, 
the right hand side of Eq. [T] gives 



4-7T a a 



2 2 

n r 



4-7T «a n 



-7T aA riQ r 3 



1 



(1 

3r c 



r/r c ) 4 
+ 3r 2 ' 



dr 



(r c + r) s 



,(7) 



which yields r c « 3.7 kpc for 2.9 x 10 56 A < 912 A 
photons per second. Note that the ceiling luminosities 
are more pronounced than for the NFW profile, due to 
the steeper power law rise at r < r c , similar to that of 
the exponential gas distribution (Fig. [SI top panel). For 
a spherical mass distribution the profile gives 



M, 



gas 



47r no 



(1 



= 47r no r c 



1 + 



- 2 ln(r 



- r/r c ) 2 
21n(r c ) 



dr 



(8) 



which again gives an physically unrealistic total gas mass, 
Mgas « 9 x 10 16 M . 

Based on the derived gas masses, we therefore con- 
clude that the exponential distribution is the most rele- 
vant to the density profile of the CNM within the host 
galaxy, although the alternative models do reproduce 
the observed critical A < 912 A luminosity for a simi- 
lar scale-length (all of which, at w 3 kpc, are close to 
that of the Milky Way). Although the e xponential dis- 
tribu tion may apply to our own galaxy (Kalbcrla et al. 
I2007T ) and, through the assumption of the canonical val- 
ues, reproduces the observed critical luminosity and gas 
mass, currently we do not have the means to map the 
density profile of the CNM at high redshift. Saying 



7 Sin ce a oc VT (|http://amdpp.phys.strath.ac.uk/tamoc/DATA/RR/[ |this, the important point is that all of these distribu- 

8 See lCurran et al.l H2UU8al) tor an inventory of the various masses tions introduces a Critical luminosity above which all of 

in a near-by active galaxy. the gas is ionised: At lower radii, when the last term 
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Fig. 4. — The rate of ionising (A < 912 A, bottom ) photons versus redshift for the sources searched in 21-cm absorption, obtained from 
the polynomial fits to the SEDs jCurran ct al. 20T2T), The symbols and histograms are as per Fig. [T] 
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n e 



-r/R 



with the top panel showing the radius of ionised gas for each value 
of R. 

in Eqs. [2j [5] and [7] are important, a power law is seen 
where log 10 r- lon oc Iog 10 J(L v /hf) dv, which is the also 
the c ase for a constant gas density profile (|Osterbrockl 
Il989t) . However, for a sufficiently high luminosity, this 
term is negligible and rj on — > oo, where the decreasing 
UV flux with distance from the AGN maintains a suffi- 
cient number photons to fully ionise the thinning gas. 

3. DISCUSSION AND CONCLUSIONS 

From the fir st high redshi f t surve y of associated 21- 
cm absorption, iCurran et al.1 (|2008b[ ) found an apparent 
lack of cool neutral gas within the hosts of radio galaxies 
and quasars. Upon an analysis of the photometry of each 
source, a strong correlation between the A = 1216 A con- 
tinuum luminosity of the AGN and the non-detection of 
21-cm was found. Specifically, that 21-cm absorption has 
never been detected in a source in which the luminosity 
exceeds L1216 ~ fO 23 W Hz -1 . Although other factors 
may contribute to the raising of the spin temperature of 
the gas, the fact that 21-cm cannot be detected above 



o 

CD 



o 

° cn 

Q. 

«: 

2 o 

CO 10 
V 



-2 °. =-*ri;-H-ti-Hr - - i . i n -; aj . j. 1 : n ;i!| — 1 1 1 11111 



C CD 



O O 

d" 
* "o 

CN J - 
i__i 01 

o 

co ,— 
"> », 

E ° 

go 
o «- 



I — I I I 1 1 III 



n = 10e 



T 1 ■ I ■ 1 1 1 1 

-r/R 



-| 1 I I I I III 



- 1 1 I I I 1 1 1| 

R = 



-I — I I I I I III' 



10000 



1000 



100 



10 



R = 1 



I I I M ill 



\ 



\ 



\ 



\ 



\ 



I hi 1 1 1| 



\ 

i i rtn ii 



\ 

I i' l H Il 



/ 



i r i t -iiii 



/" . - 





1 



10 1 00 1 000 

Galactocentric radius, 



10° 



10" 
r [pc] 

Fig. 6. — As Fig. [5]but for a gas temperature of 10 000 K (oa = 
4.19 X 10 -13 cm 3 sec -1 ) and a spherical gas distribution. 

the ground state, in conjunction with the lack of detec- 
tions above this critical luminosity (significant at 4.76cr), 
strongly suggests that excitation byA<1216A photons 
is the dominant cause of the dearth of 21-cm in optically 
bright radio sources. 

Here we demonstrate that this critical luminosity is 
also applicable to ionising (A < 912 A) photons, show- 
ing that associated 21-cm is not detected for any source 
where L912 > 10 23 W Hz -1 or, more precisely, when there 
are Z 2.9 x 10 56 ionising photons sec -1 . Applying this 
photoionsation rate, together with various gas density 
distribution models to the equation of photoionsation 
equilibrium, from canonical values for the gas density 
(no = 10 cm -3 ) and the recombination rate coefficient 
(a A = 1-27 x 10~ 12 cm 3 sec" 1 ): 

• We obtain the observed critical photon rate for a 
scale-length of w 3 kpc for all of the tested profiles 
(exponential, NFW and isothermal sphere). This 
scale-length is the same as that for the Hi in the 
Milky Way, thus suggesting that the observed criti- 
cal value is just sufficient to ionise all of the neutral 
gas within a large spiral galaxy. 

• This scale-length gives: 
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Fig. 7.— As Fig. [5] (no = 10 cm" 3 and T = 2000 K) but for an 
NFW profile. 

— For an exponential distribution within a disk, 
a total gas mass of M gas = 7.5 x 10 9 M , typ- 
ical of that found from 21-cm emission studies 
of low rcdshift galaxies. 

— For the NFW and isothermal sphere distribu- 
tions, a total gas mass which exceeds the total 
expected dynamical mass of the galaxy. 

This leads us to conclude that the exponential profile is 
the more applicable to the distribution of the CNM (see 
iKalberla fc Kerdl2009D . although all of the models give 
a critical UV luminosity. That is, for a gas profile in 
which the density decreases with distance from the nu- 
cleus, the Stromgren sphere has an infinite radius for a 
finite luminosity. This suggests that a balance is main- 
tained between the decreasing number of photons and 
number of particles with increasing distance from the 
ionising source. 

For the sources under consideration here, the criti- 
cal photon rate (where L912 ~ £1216 ~ 10 23 W Hz -1 ) 
is consistent with the dearth of 21-cm detections i n 
all searched high rcdshift sources (jCurran et al.ll2008bD . 
A "proximity effect" for highly i onised Lyman-q forest 
clouds has previou sly been noted dWevmann et alj|1981t 
iBaitlik et al.l 119881 ). where in these intervening systems 
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Fig. 8.— As Fig. [7j (ra = 10 cm" 3 and T = 2000 K) but for a 
halo profile with an isothermal sphere density distribution. 

the high ionising flux from the QSO is believed to be re- 
sponsible for the decrease in the number density of the 
Lyman-a lines as the rcdshift of the absorbing galaxy 
approaches that of the QSO (z a b s — > zqso)0 

However, until our high redshift survey of radio galax- 
ies and quasars, no such effect was known for the 21-cm 
transition! 1 "! Furthermore, the 21-cm effect is striking in 
that, rather than a gradual decrease in associated 21-cm 
absorption with increasing ultra-violet luminosity, there 
is an abrupt cut-off in the 21-cm detection rate at a single 
critical luminosity. We show here that, for a typical spec- 
tral energy distribution, this is the luminosity required 
to completely ionise the large-scale distribution of atomic 
gas in a large spiral galaxvFH 

9 IBah call & Ekcrs (1969) also show that both the 21-cm and 
the Lyman-a flux can contribute to higher spin temperatures at 
absorber-quasar separations of less than a few tens of kpc. 

10 From a study of ab sorber clustering around QSOs in the SDSS 
DR3, IWild et alt 020081 ) suggest that the QSO destroys the Mg II 
clouds out to beyond 800 kpc. Mgll has an ionisation potential 
of 15.04 eV, close that of Hi (13.60 eV) and so a similar critical 
A < 827 A luminosity could perhaps account for this. 

11 Although not searched in 21-cm, we derive an ionising pho- 
ton rate of 3.6 X 10 57 sec" 1 (L912 = 2.1 X 10 24 W Hz" 1 ) for 
PKS 0424-131. Since this is an order of magnitude higher than 
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There is the possibility that the non-detection of 21-cm 
absorption in the high-luminosity sample is simply due to 
orientation effects - as mentioned previously (SectfTJ, the 
majority (17 of 19) of the Luv ~ 10 23 W Hz -1 sources 
are type-1 AGN. In the plane of the torus there may exist 
large columns of neutral gas, shielded from the radiation 
(|Hennawi fc Procha ska 2007). The orientation of these 
objects could explain their high UV luminosities, how- 
ever there are also many type-1 objects below Luv ~ 10 23 
W Hz -1 , which exhibit a 50% probability of detecting 
21-cm absorption. This suggests that the absorbing gas 
is located in the large-scale disk, which must be ran- 
doml y oriented with re s pect to the torus (iSchmitt et all 
T997t INaear fc Wilson! [19991: [C urran fc Whitind I20TqT 
Lawrence fc Elvidl2O10l ). 

Therefore, if the non-detection of absorption in the 
Luv ~ 10 23 W Hz -1 sources is due to line-of-sight effects, 
these differ from the lower luminosity AGN in that, un- 
like these, the large-scale absorbing disk is always aligned 
with the obscuring torus. It is possible that for these, the 
high UV luminosities are indeed the consequence of a di- 
rect view to the AGN unobscured by the large-scale disk. 
This would also explain the lower luminosity type-1 ob- 
jects in which 21-cm is detected, in that the absorbing 
gas attenuates the radiation. 

Using the orientation as the sole explanation for the 
differences in absorption rates can not, however, explain 
why there is a critical value in the UV luminosities rather 
than a gradual decrease in absorption rates as the lu- 
minosity increases. The model presented here naturally 
yields a critical luminosity, and we have shown that this 
value is sufficient to ionise all of the CNM in a large spiral 
galaxy. 

Our model is supported by the observation of greatly- 
decreased 2 50^m emission f rom AGN above a critical 
luminosity (|Page et alll2012[) . The emission is not re- 
stricted to the line-of-sight to the nucleus, indicating that 
the entire galaxy is affected by the presence of a luminous 
AGN. 



Finding rcdshiftcd 21-cm absorption is a major science 
goal of the Square Kilometre Array (SKA). The results 
found here suggest that the lack of cool neutral gas is not 
due to a sensitivity issue with current radio telescopes, 
but a real effect caused by the presence of a quasar in 
a galaxy of gas. That is, the known high redshift radio 
galaxies and quasars are probably devoid of a large-scale 
distribution of neutral gas\^\ This neutral gas provides 
the fuel for star formation, and our result suggests that 
the AGN could therefore suppress star formation in the 
large-scale disk. Note also that, although the 250 /im 
fluxes are also subject to a sensitivity limit, the deficit of 
250 fim above a critical luminosi ty is further evid ence of 
a suppression of star formation (Page et al. 2012). Here 
we show that this is not a sensitivity issue, but that the 
neutral gas, and most likely, appreciable star formation 
activity is simply not present. 

Therefore, even the SKA will be unlikely to find this 
cool gas in the objects currently known. Where it will 
excel, however, is in blind surveys of radio sources from 
which the visible lig ht is too faint to b e detected by 
optical instruments (ICurran etaLl [2009T) . Although 
unseen, these sources must exist in order to have had 
star formation within the host galaxies of early AGN. 
As such, the traditional optical selection of targets must 
be abandoned in order to find the missing star-forming 
material within high redshift radio sources. 
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vice. The Centre for All-sky Astrophysics is an Aus- 
tralian Research Council Centre of Excellence, funded 
by grant CE1 10001020. 



REFERENCES 



Allison, J. R., Curran, S. J., Emonts, B. H. C, et al. 2012, 

MNRAS, 423, 2601 
Athreya, R. M., Kapahi, V. K., McCarthy, P. J., & van Breugel, 

W. 1998, A&A, 329, 809 

the critical rate, it could account for "the mysterious absence of 
neutral hydrogen", as traced by Lyman -c t emission, c l ose to this 
source flf^ancis fc Bland-Haw thorn 200D. IBruns et"aTI H2012I) also 
suggest that the UV emission could be responsible, thus suppress- 
ing the star formation. 

12 Simulations suggest that, while the Hi content in galaxies 
is similar to present day values, the molecular component (H2) is 
larger at z Z 1 (Obrcschkow fc Rawlings 2009). For gas densities 
typical of molecular clouds, no ~ 1000 cm~ 3 , the critical ionising 
photon flux of 2.9 X 10 56 sec -1 gives a scale-length of a 100 pc at 
2000 K, or ss 50 pc at a more realistic 20 K. Conversely, luminosities 
of L912 w 2xl0 27 WHz" 1 (3xl0 60 photons sec" 1 ) and4xl0 28 W 
Hz" 1 (5 X 10 61 photons sec" 1 ) are required to fully ionise a CNM 
with no ~ 1000 cm -3 with a scale-length of 3 kpc at 2000 and 
20 K, respectively. If the deficit of neutral gas at high redshift was 
due to the gas being mostly molecular, there would be little reason 
for a lack of 21-cm to be correlated with the A = 912 A luminosity 
(although a lack of H2 could be). In any case, the Hi deficit is 
observed at all redshifts, as well as there being no observational 
evidence for a high m olecular gas content w ithin the hosts of radio 
galaxies and quasars ICurran et al.ll2011bl) . 



Bahcall, J. N. & Ekers, R. D. 1969, ApJ, 157, 1055 
Bajtlik, S., Duncan, R. C, & Ostriker, J. P. 1988, ApJ, 327, 570 
Begeman, K. C, Broeils, A. H., & Sanders, R. H. 1991, MNRAS, 
249, 523 

Bruns, Jr., L. R., Wyithe, J. S. B., Bland-Hawthorn, J., & 

Dijkstra, M. 2012, MNRAS, 2543, 421 
Catinclla, B., Haynes, M. P., Giovanelli, R., Gardner, J. P., & 

Connolly, A. J. 2008, ApJ, 685, L13 
Curran, S. J. 2010, MNRAS, 402, 2657 
— . 2012, ApJ, 748, L18 

Curran, S. J., Koribalski, B. S., & Bains, I. 2008a, MNRAS, 389, 
63 

Curran, S. J., Tzanavaris, P., Darling, J. K., et al. 2010, MNRAS, 
402, 35 

Curran, S. J. & Webb, J. K. 2006, MNRAS, 371, 356 
Curran, S. J., Webb, J. K., Murphy, M. T., et al. 2002, PASA, 19, 
455 

Curran, S. J., Whiting, M. T., Murphy, , M. T., et al. 2011a, 

MNRAS, 413, 1165 
Curran, S. J. & Whiting, M. T. 2010, ApJ, 712, 303 
Curran, S. J., Whiting, M. T., Combes, F., et al. 2011b, MNRAS, 

416, 2143 

Curran, S. J., Whitin g, M. T., Sadler, E. M., & Bignell, C. 2012, 

MNRAS, in press darXiv:1210.1886t 
Curran, S. J., Whiting, M. T., Tanna, A., Bignell, C, & Webb, 

J. K. 2011c, MNRAS, 413, L86 



Complete ionisation of the neutral gas in high rcdshift radio galaxies and quasars 



13 



Curran, S. J., Whiting, M. T., & Webb, J. K. 2009, Proceedings 

of Science, 89, Chap. 11 
Curran, S. J., Whiting, M. T., Webb, J. K., & Athrcya, A. 2011d, 

MNRAS, 414, L26 
Curran, S. J., Whiting, M. T., Wiklind, T., et al. 2008b, MNRAS, 

391, 765 

de Blok, W. J. C, Walter, F., Brinks, E., et al. 2008, AJ, 136, 
2648 

Draine, B. T. 2011, Physics of the Interstellar and Intergalactic 

Medium (Princeton University Press) 
Field, G. B. 1959, ApJ, 129, 536 

Francis, P. J. & Bland-Hawthorn, J. 2004, MNRAS, 353, 301 

Copal-Krishna & Wiita, P. J. 2000, ApJ, 529, 189 

Grasha, K. & Darling, J. 2011, in American Astronomical Society 

Meeting Abstracts, Vol. 43, 345.02 
Gupta, N., Salter, C. J., Saikia, D. J., Ghosh, T., & Jcyakumar, 

S. 2006, MNRAS, 373, 972 
Haiman, Z. & Rees, M. J. 2001, ApJ, 556, 87 
Hennawi, J. F. & Prochaska, J. X. 2007, ApJ, 655, 735 
Jaffe, W. 1983, MNRAS, 202, 995 

Kalberla, P. M. W., Dedes, L., Kerp, J., & Haud, U. 2007, A&A, 
469, 511 

Kalberla, P. M. W. & Kerp, J. 2009, Ann. Rev. Astr. Ap., 47, 27 
Koribalski, B. S., Staveley-Smith, L., Kilborn, V. A., et al. 2004, 
AJ, 128, 16 

Lawrence, A. & Elvis, M. 2010, ApJ, 714, 561 

Martin, C, Barlow, T., Barnhart, W., et al. 2003, in Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 4854, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, ed. J. C. Blades & 
O. H. W. Siegmund, 336-350 

Morganti, R., Oosterloo, T. A., Tadhunter, C. N., et al. 2001, 
MNRAS, 323, 331 

Nagar, N. M. & Wilson, A. S. 1999, ApJ, 516, 97 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 
563 

Noterdaeme, P., Petitjean, P., Ledoux, C, & Srianand, R. 2009, 
A&A, 505, 1087 



Obreschkow, D. & Rawlings, S. 2009, ApJ, 696, L129 

Oh, S.-H., de Blok, W. J. C, Brinks, E., Walter, F., & Kennicutt, 

Jr., R. C. 2011, AJ, 141, 193 
Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and 

Active Galactic Nuclei (Mill Valley, California: University 

Science Books) 

Osterbrock, D. E. & Ferland, G. J. 2006, Astrophysics of gaseous 
nebulae and active galactic nuclei (Sausalito, California: 
University Science Books) 

Page, M. J., Symeonidis, M., Vieira, J. D., et al. 2012, Nature, 
485, 213 

Peroux, C, Storric-Lombardi, L. J., McMahon, R. G., Irwin, M., 

& Hook, I. M. 2001, AJ, 121, 1799 
Pihlstrom, Y. M., Conway, J. E., & Vermeulen, R. C. 2003, A&A, 

404, 871 

Purcell, E. M. & Field, G. B. 1956, ApJ, 124, 542 
Schlcgel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 
525 

Schmitt, H. R., Kinney, A. L., Storchi-Bcrgmann, T., & 

Antonucci, R. 1997, ApJ, 477, 623 
Scott, J. E., Kriss, G. A., Brotherton, M., et al. 2004, ApJ, 615, 

135 

Shull, J. M., Stevans, M., & Danforth, C. W. 2012, ApJ, 752, 162 
Srianand, R. , Gupta, N., Petitjean, P., Noterdaeme, P., & 

Ledoux, C. 2010, MNRAS, 1888 
Taylor, G. L., Dunlop, J. S., Hughes, D. H., & Robson, E. I. 1996, 

MNRAS, 283, 930 
Telfer, R. C, Zheng, W., Kriss, G. A., & Davidson, A. F. 2002, 

ApJ, 565, 773 

Uson, J. M., Bagri, D. S., & Cornwell, T. J. 1991, PhRvL, 67, 
3328 

Weymann, R. J., Carswell, R. F., & Smith, M. G. 1981, ARA&A, 
19, 41 

Wild, V., Kauffmann, G., White, S., et al. 2008, MNRAS, 388, 
227 

Wolfe, A. M. & Burbidge, G. R. 1975, ApJ, 200, 548 



